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Abstract. Persistent homology was shown by Zomorodian and Carlsson | 35| to be homology of graded 
chain complexes with coefficients in the graded ring k[t]. As such, the behavior of persistence modules — 
graded modules over k[t] — is an important part in the analysis and computation of persistent homology. 

fn this paper we present a number of facts about persistence modules; ranging from the well-known 
but under-utilized to the reconstruction of techniques to work in a purely algebraic approach to persistent 
homology. In particular, the results we present give concrete algorithms to compute the persistent homology 
of a simplicial complex with torsion in the chain complex. 
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1. Introduction 

The ideas of topological persistence [20] and persistent homology [35[ have had a fundamental impact on 
computational geometry and the newly spawned field of applied topology. The method cleverly modified 
tools from algebraic topology to make them resistant to noise and applicable to input from applications. 
It has found many uses, and produces globally descriptive results often completely inaccessible with other 
methods. Theoretical work in the field has served to successively strengthen the weight and applicability of 
the methods. For surveys or books on the subject, we recommend: [fj, 0, HH, Efijl. 

Since the initial publication Edelsbrunner, Letscher, and Zomorodian [ioj ]. persistent homology has spawned 



many results in applied and computational topology. It has been applied to manifold learning [3CJ, |31| . 
bioinformati cs [lil liHj. computational chemistry 27 1, medical @,[2^|, trace data [II], among others (see also 



0, H E > 12, 32|) . See [l9[ for a complete survey of the development the field 

The introduction of algebraic insights into the field has already been productive: Zomorodian and Carlsson 
[35| identified persistence modules as graded modules over the polynomial ring \k[t] over a coefficient field k, 
and used this insight to produce advances in the computation of persistent homology. Even later, Carlsson, 
de Silva, and Morozov [9( recognize that the tame representation theory of A„-quivers allow an analogous 
construction producing zig-zag persistence. 

Building on this backdrop, we introduce techniques from computational commutative algebra into the study 
of persistent homology - maintaining explicit presentations of persistence modules enables us to give purely 
algebraic approaches to kernels, cokernels and images of morphisms between persistence modules - providing 



an alternative to the approach by Cohen- Steiner, Edelsbrunner, Harer, and Morozov 15j. In addition to 
this, it provides algorithms for a variety of other algebraic constructions including algorithms to compute 
persistent homology of chain complexes with torsion present, and on that route an algebraic approach to 
relative persistent homology and cohomology. 



1 .1 . Prior work. At the core of the study of topological persistence lies two papers: Edelsbrunner, Letscher, 
and Zomorodian [20| who introduced the idea of persistent homology in the first place, working exclusively 
with coefficients in Z/2Z and giving an ad hoc algorithm for computing persistence barcodes. 

This was followed by Zomorodian and Carlsson (35|, who identified persistence modules as, essentially, graded 
modules over the graded polynomial ring k[t], and described based on this identification how algorithms follow 
where k may be an arbitrary field. 



The algebraic formalization approach, combined with results from Gabriel [22] produced fundamental results 
by Carlsson, de Silva, and Morozov [9(, identifying quiver algebras as interesting coefficient rings, producing 
a theory of zig-zag persistent homology that allows the study of diagrams of topological spaces that no longer 
form strict nitrations. As such, these form one way to approach the question of how to handle torsion, or 
vanishing simplices, in a persistent homology approach. 

The usefulness of persistent homology relies fundamentally on stability results - for persistent homology 
generated from sublevel filtrations of topological spaces, a small change in the filtration function generates 
a quantifiably small change in the resulting persistence modules. These types of results were introduced to 



the field by Cohen- Steiner, Edelsbrunner, and Harer [14{ who require the underlying space to be triangulable 



and the persistent homology modules to be tame and degree-wise finite dimensional. Their results were 
significantly improved by Chazal, Cohen-Steiner, Glisse, Guibas, and Oudot (To| . who formalize persistent 
homology as R-indexed diagrams of vector spaces, and define an interleaving distance between such diagrams 
to provide stability results that no longer require continuous filtration functions, triangulable spaces, or tame 
persistent homology modules. These stability results culminate in recent work by Bubenik and Scott [Bj], 
who study the category of R-indexed diagrams in arbitrary abelian categories, and are able to prove a 
range of stability results on sublevel filtrations using arbitrary filtration functions for arbitrary functors from 
topological spaces to some abelian category. 
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Algebraic operations in persistent homology have also been considered before. The kernel, cokernel and image 
constructions we introduce here are alg ebraic reformulations of the techniques introduced in Cohen-Steiner, 
Edelsbrunner, Harer, and Morozov [15j . 

Torsion in chain complexes is tightly linked to relative homology, a tool that has shown up numerous times 
in the field. Cohen-Steiner, Edelsbrunner, and Harer provide a framework for persistence in which by 
including relative homology groups of the entire space relative to a superset filtration, all barcode intervals arc 
finite in length. This approach is refined by Carlsson, de Silva, and Morozov [7] , who are able to fit extended 
persistence into a Mayer- Vietoris pyramid, connecting various constellations of relative persistent homologies 
into a large and expressive system. Bendich, Cohen-Steiner, Edelsbrunner, Harer, and Morozov Bendich, 
Wang, and Mukherjee [3], and Wang, Summa, Pascucci, and Vejdemo- Johansson [34[ use persistent versions 
of local homology - homology relative the complement of a small neighborhood of a point - to discover local 
behavior close to that point. 



1.2. Our contributions. In this paper, we introduce techniques from computational commutative algebra 
to the study of persistence modules. In Section we review the correspondence between persistence mod- 
ules and graded Ik[i]-modules. Section 13.11 reviews fundamental concepts from computational commutative 
algebra. Of importance to the rest of the paper is the adaptation of Smith normal form in Section 13.1.31 
which gives a graded version of the Smith normal form and its computation that we have been unable to 
find in the literature. Several of the results in Section 13.11 are regularly discussed either for simpler rings 
(fields, integers) or vastly more complex rings (algebraic geometry) - we discuss the special case of graded 
lk[t]-modules and their applications to persistence. In particular, the maximum possible free dimension in the 
category of graded lk[i]-modules is discussed with its implications for persistence modules. The goal of this 
paper is twofold: to expand on the algebraic interpretation of persistence and provide a uniform framework 
for algebraic constructions using persistence modules 

In Section I3.2[ we start elaborating on standard techniques from computational algebra for representing 
finitely presented modules, and draw from these to give specific constructions for persistence modules. Sec- 
tion discusses how a presentation map from relations to generators captures the behavior of a persistence 
module, and allows a graded Smith normal form computation to recover a barcode from an arbitrary pre- 
sentation. The sections 13.51 13.61 13.71 13.81 13.101 and 13.111 are dedicated to purely algebraic constructions of 
kernels, cokernels, images, pushouts, and pullbacks of persistence module maps. 

In Section IU we discuss applications of the techniques and algorithms that we have described. In Section |4T| 
we show how the nested module presentation of a quotient module together with kernel and cokernel algo- 
rithms allows us to compute persistent homology of a chain complex with torsion, giving a purely algebraic 
approach to relative homology. 

In Section [4.21 we discuss the relations between the work in this paper, and forthcoming work on computing 
spectral sequences of persistence modules. 

In Section [4.3l we demonstrate how the fundamental viewpoint of graded modules over a graded ring provides 
us a variation of the persistence algorithm that takes an unsorted stream of simplices as input, and changes 
the output on the fly should new simplices indicate that the previous inferences were inaccurate. 



2. Persistent (Co-)Homology 

In persistent homology and cohomology, the basic object of study is a filtered simplicial complex, with the 
filtration often induced as a sublevel- or superlevel-filtration of some (real- valued) function on a space. Since 
the filtration can be considered a special kind of representation in simplicial complexes of the total order of 
the function domain, where all the endo-maps induced in the representation by the total order represented 
are injective, it is clear that because of functoriality, we can take point-wise homology (with coefficients in 
Ik) on the representation and thus get a functor 

SpCpx Domain -> ModJ? omain 
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taking each filtered simplicial complex to a directed system of Ik-modules, in which the map H*(X e ) — > 
H*(X e i) is induced by the inclusion X e X e > for e < e'. 

For finite simplicial complexes, it is clear that at most finitely many critical points Xi can exist such that 
for each Xi, and for all sufficiently small e, the map X Xi _ E <^-» X Xi+e are not isomorphisms. Hence, we can 
change the entire representation from the original range to a representation over a finite order 

—oo < (xq + x\)/2 < [x\ + xa)/2 < • • • < (x n _i + x n )/2 < oo 

The same observation holds for a countable simplicial complex, which corresponds to a countable order 
bounded by the top and bottom ±00. 

Assuming that the order is left bounded, or in other words, that for any point along the order there are only 
finitely many populated "compartments" less than the current one, there are clear isomorphisms between 
chain complexes for these representations and graded modules over the ring k[t]. Indeed, given left bound- 
edness, we can write for the center points stated above of the compartments ao = —00, a\ = (xq +&i)/2, . . . 

enumerating the compartments as living above ao, ai, We introduce a N-grading on the chain complex 

by making an element that lives over au have degree k. Finally, we introduce a lk[f]-action by letting 

(1) r ■ c be given by the Ik-module structure on the chain complex C„X ah for r £ Ik and c £ C*X ak 

(2) t ■ c be given by the induced map C*X aic — > C^X ak+1 

This extends by lk[i]-linearity to a 0<[t]-module structure for the entire representation C,X,. 

At the core of the observations by Zomorodian and Carlsson [35| , and also by Carlsson, de Silva, and Morozov 
Q is that the computation of persistent homology is really the computation of homology within the category 
of modules over the appropriate quiver algebra - D<[£] or even lk[i]/i n if we have classical persistence, and 
kQ for a quiver Q of type A n for zig zag persistence. This observation underlies our approach to algebraic 
methods in persistent homology. 

Indeed, once we can identify C*X» with a k[i]-module by the methods above, the remainder of a persis- 
tent homology computation follows immediately. The persistence algorithm as described by Zomorodian 
and Carlsson [35| is a Gaussian elimination algorithm as applied to graded (k[i]-modules, and the entire 
computation follows from this recognition. 

3. Overview of Commutative Algebra 

In this section, we will recall important and relevant facts from commutative algebra; for some of these 
results we will be using the internals of the proofs in later developments. The entirety of Section [3. II can be 
skipped if you are already comfortably familiar with algebra. 

Commutative algebra is an area of study which deals with commutative rings and their associated objects. 
With a field, the lessons learned from commutative algebra reduce to classical linear algebra - and so, for 
much of what we will be dealing with, the commutative algebra can well be considered to be a generalization 
of familiar techniques from linear algebra into areas where not all the strength of linear algebra is available. 

In particular, the construction of a basis for an arbitrary finite dimensional vector space is a fundamental 
building block in many algorithms. It, however, breaks down if the coefficients do not lie in a field. The core 
step is where the leading coefficient in some potential basis vector gets reduced to I by dividing all coefficients 
with this leading coefficient. In general commutative rings, there is no guarantee for having multiplicative 
inverses, and therefore, this step can fail. 

Instead, in commutative algebra, for a vector space like object - a module - to have a basis is a very particular 
feature of the module, and is enough to warrant a separate name for that kind of module: a free module. 
As we shall see later on, a lot of the power comes from phrasing our questions and our objects in ways that 
retain free modules in the descriptions and solutions as much as possible. 

3.1. Ring Theory. A ring is some set of elements together with two binary operations: addition (denoted 
by +) and multiplication (denoted by *, • or juxtaposition). These are made to follow reasonable axioms: 
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addition should make the entire ring into an abelian group. In particular, addition is commutative, associa- 
tive, and for each element a, it has an additive inverse —a, as well as an additive identity clement in the 
ring. 

Multiplication also follows axioms: associativity, and commutativity (hence commutative algebra), and we 
shall also always require a multiplicative identity element (denoted by 1). 

Finally, the two operations are tied together by distributivity: (a + b) ■ c = a ■ c + b ■ c. 

The notable difference from the well-known definition of a field is the lack a requirement of multiplicative 
inverses: indeed, if each non-zero element a has a multiplicative inverse a -1 such that a ■ a^ 1 = 1, then the 
commutative ring is in fact a field. 

A "vector space" M over a commutative ring R is called a module (or more accurately, a vector space 
is a module over a field). More precisely, a module is an abelian group with a specific binary operation 
• : R x M — > M which obeys the following axioms: 

• • m = 0, 1 • m = 1, 

• r • (s • m) — (r ■ s) ■ m; (r + s) ■ m — r ■ m + s ■ m, 

• r ■ (m + n) = r ■ m + r ■ n for r, s G R and m,n G M. 

We now state a few technical definitions followed by some informal comments. 

A submodule N < M is a subgroup N C M such that for any r G R and any n G N, r ■ n G N. This is 
analogous to subspaces of vector spaces. The submodule is trivial if it only contains the element, and it is 
proper if it is a strict subset of M. Given a submodule N < M, the quotient group M/N has the structure 
of an i?-module by defining r • [m] — [r ■ m]. We call the resulting module the quotient module M/N. A 
presentation of a module M is a presentation of M as a quotient module F/K. The presentation is free 
if both F and K are free modules. The module M is finitely generated if a presentation exists with F 
free with a finite basis. The module M is finitely presented if both F and K are free modules with finite 
bases. 

A submodule N < R oi the ring itself is called an ideal. 

A ring is graded over a monoid (G, +) if the ring decomposes as a direct sum R = © geG R g where each 
R g is a subgroup of the additive group of R, and for elements r G R g and s G Rh, the product is in R g+ h- A 
module M over a graded ring R is graded if it decomposes into a direct sum of subgroups M = © ggG M g 
such that the scalar product obeys the grading: r G R g and m G Mh implies r ■ m G M g+ h- 

A ring is noetherian if any ascending chain of ideals Iq C I x C . . . eventually stabilizes. In other words, R 
is noetherian if for every such chain of ideals there is some integer N such that for all n > N, I n = I n +i- 
Very many rings we know are noetherian, in particular Z, all fields, all principal ideal domains, and all 
polynomial rings over noetherian rings. Any ideal in a noetherian ring is finitely generated, and the partially 
ordered (under inclusion) set of any non-empty set of ideals of a noetherian ring has a maximal element. 

Any ring has at least one proper ideal: the set {0} is always an ideal, usually denoted by with some abuse 
of notation. If is the only proper ideal in a ring i?, then that ring is a field. Indeed, suppose Then, 
since is the only proper ideal, the principal ideal (r) has to include the entire ring. In particular, I G (r), 
and thus there is some s such that rs = 1. Hence, r has a multiplicative inverse. 

Theorem 1 (Decomposition of modules over PIDs). Suppose R is a principal ideal domain. Then for every 
finitely generated module M over R, there is an integer n and a sequence of elements ri,...,r m such that 




The proof of this is standard, and relies on the observation that matrices over a principal ideal domain admit 
a Smith normal form. 
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Figure 1. Persistent homology example by Edelsbrunner, Letscher, and Zomorodian [201 ] . 
The arrow denotes the direction of time in the filtration with the graded complex shown at 
the far right. 



3.1.1. Euclidean Domains. The following results and their proofs are well-known facts of commutative and 
computational algebra. We shall, however, reproduce some of them in order to be able to refer to their 
particular structure in this paper. 

Definition 2. Suppose R is a ring. A Euclidean function on R is a function R — > N such that if a,b £ R 
and b ^ 0, then there are q,r £ R such that a = b ■ q + r and either r = or deg r < deg b. 

An integral domain R that supports at least one Euclidean function is called a Euclidean domain. 
Theorem 3. In a Euclidean domain, any two elements have a GCD. 

Proof. The extended Euclidean algorithm relies precisely on the existence of a Euclidean function, and 
generates, for elements r, s elements a, b, g such that g = gcd(r, s) and a ■ r + b ■ s = g. □ 

Theorem 4. All Euclidean domains are PIDs. 

Proof. Suppose I C R is an ideal in a Euclidean domain. By well-order of N, there is some minimal value in 
the set {degr|r £ /}, say n. Again, by the well-ordering of N, this value is attained - so there is some b £ I 
such that deg b = n. 

We claim that / is generated by b. Indeed, suppose that a £ / is some other element of /. Since deg a > deg& 
by the minimality of deg b, there is some q,r such that a — b ■ q + r and either r = or degr < deg b. Since 
r = a — b ■ q, and both a, b £ I, it follows that r £ I. 

Hence, if r ^ 0, then degr < deg& contradicting the minimality of degree in our choice of b. Hence, r = 
and thus a = b ■ q, so a £ (b). 

The result follows immediately. □ 

3.1.2. Connection with Persistence. Zomorodian and Carlsson [35j j identify the persistence modules with 
graded modules over k[t\. We remind the reader of the intuition of this representation. Recall that persistent 
homology computes the homology over a nested sequence of spaces called a filtration: 

= X C Xi C X 2 C . . . C X w _i C Xat = X 

The grading comes from viewing each of the spaces in the sequence as a slice and stacking the slices up as 
in Figure [TJ 

This entire structure is then encoded in the final complex (top slice), with each simplex annotated by a 
number indicating how many slices ago it first appeared. This is precisely the information encoded in the 
k[t] module structure. 
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Furthermore, the free and torsion parts of the decomposition in Theorem Q] correspond to infinite and finite 
intervals in the persistent barcodes respectively. With these results in hand, we can state the following 
results. 

Theorem 5. The graded ring k[t] is a Euclidean domain. 

Proof. As a Euclidean function, take deg^a^' = max{i|ai ^ 0}. We notice that this coincides with the 
classical notion of degree of a polynomial, and with the grading of k[t]. □ 

Theorem 6. Any finitely generated submodule AI of a free module F is itself free. 

Proof. Suppose M C F is generated by elements mi, . . . ,mk- Suppose that M is not free. Then there is 
some Ik [t] -linear relation 



between the generators. If e t > whenever ^ 0, then the entire expression is a multiple of t. Since F is 
free, if t ■ m = then m = for all elements m £ F. In particular, this holds for t ■ a,it ei ~ 1 m i . Thus, we 
reduce the exponents in this expression by some integer d until for some i, ^ and e,; = 0. Thus, this 
relation is equivalent to the relation 



But this implies that rrii is in the span of mi, . . . , Wj, . . . , m/j, where we retain the convention that the rhi 
means that we exclude mi from the list. Thus, the same module M is also generated by the remaining 
generators. 

Since there are finitely many generators, this process of elimination eventually terminates. When it termi- 
nates, it must be because there are no more non-zero k[<]-linear relations to be found. 

Hence, if M is a finitely generated submodule of F, it must be free. □ 

Theorem 7. Any finite presentation K—^G—>M—>0ofa graded module M over k[i] is a free presentation 
-> K ->■ G -> M 0. 

Proof. Suppose G is a free module of generators of M, and G A M — > is the map of basis element to 
generator. K, the module of relations of the presentation is by force the free module of generators of kerp. 

Claim: The corresponding map K A G is an injective map. 

Indeed, suppose i is not injective as a degree map of graded k[t]-modules. Then there is some k ^ such 
that ik = 0. If k is not homogeneous, then k — kd, and there is some non-zero homogeneous kd such 
that ikd = 0. 

Consider the map restricted to the degree d part of all modules involved: Kd Gd Mj, — > 0. Since 
everything is graded, this restriction of a presentation of M is a finite presentation of Md as a Ik-module, 
or in other words as a vector space. Since all vector spaces are free modules, it follows that kerp^ is a free 
subspace of Gd, and thus WLOG, ker pd = Kd with id the inclusion map of kerp^ C Gd- Therefore, id is an 
injective map, and thus has a trivial kernel. It follows that idhd — implies kd = 0. But this contradicts the 
assumption above. Hence i is injective, and the result follows. □ 

It is worthwhile to demonstrate exactly how this representation works. Chains in a graded 0<[£]-module of 
chains will have an expression that depends both on the birth times of their basis elements and on the point 
in time that the chain itself is considered. Consider the chain ab + ac + cd in Figure [TJ existing at the 5th 
time step of the diagram. To find a [k[i]-linear description of this chain, we need to promote each of the three 
basis elements, ab, ac and cd, into the 5th time step. Since the edge ab shows up in the second time step, 
this takes a coefficient of t 3 . Similarly, ac has a coefficient of t and cd a coefficient of t 2 . All in all, the chain 
is t ■ ac + t 2 ■ cd + t 3 ■ ab. 
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This is also relevant for combining elements of graded k[i]-modules: each (homogeneous) element has a 
degree it lives in — and to sum different elements, ensuring a homogeneous result, the elements need to live 
in the same degree; the cycle basis from the situation in Figure [T] has as natural choices of basis elements 
the two cycles z = ad — cd — t ■ ah — t ■ be and w = ac — t 2 ■ ab — t 2 ■ be. To express the other small cycle 
at time step 4, we will need to form a linear combination of the two cycles at hand at the time step we are 
interested in. Hence, this other triangular cycle has the expression w — t ■ z = ac — t ■ ad — t ■ cd. 



3.1.3. Graded Smith normal form. The Smith normal form of a function M — > N will compute a simulta- 
neous basis change for M and N such that in the new basis, each basis element for M pairs with a basis 
element in N; the map takes any basis element to a scalar multiple of its paired element. 

While developed for matrices over Z or N, the constructions all hold over any PID - the Smith normal form 
is one of the standard ways to prove the decomposition Theorem [1] We shall describe how it is usually 
computed, with special attention to the shortcuts and differences that the case of graded modules over k[t] 
introduce. We have been unable to find such an adaptation in the literature. 

The usual algorithm allows the user to permute both rows and columns of the matrix under treatment; since 
the grading is important to us, we shall instead mark rows and columns as finished - rather than to permute 
them. 

We represent a map M — > N by a matrix F = (fij) such that f(mi) = fijTij, where nit and rij are basis 
elements in M and N respectively. 

A basis change that replaces rij by rij — rnj> has the effect of adding r times the jth row to the j'th row. 
Similarly, a basis change replacing m,i by — rm? adds r times the ith column to the z'th column. We will 
be performing these operations from left to right, from bottom to top, in a matrix that has been pre-sorted 
to keep both bases of M and N in ascending degree order. 

The resulting algorithm is essentially forced by the requirement that we can only add basis multiples of 
compatible orders: if |mj| < |tjv|, we cannot influence to$ using 7?v - we can only influence m<> using m^. 
Hence, the entries in the matrix can only flow upwards in their influence; never downwards. 

Algorithm: 

(1) While there are untreated rows or columns: 

pick the lowest degree entry of the lowest degree block of untreated columns. 

(2) By appropriate basis changes in M, clear out all entries of column i. Entries are cleared out upwards 
in the matrix. 

(3) By appropriate basis changes in N, clear out all entries of row j. Entries are cleared out leftwards 
in the matrix. 

(4) Mark row and column as finished. 

(5) Goto step HI 

For extra clarity, we include two examples here. Consider the sequence in Figure [T] We can represent the 
chain complex by a graded k[i]-module with chain basis a,b,c,d,ab,bc,ad,cd,ac,abc,acd and boundaries 
tb — ta,c — tb, td — tc, td — t 2 a, t 2 c — t 3 a, t 3 ab + t 3 bc — tac, t 3 cd — t 3 ad + t 2 ac. Computing a kernel of the 
boundary map, we get a basis for the cycle module given by: 



z 2 =b 



Z\ = a 



z 4 = d 



z 3 = c 



z 5 — tab + tbc + cd — ad 
zq = t 2 ab + t 2 bc — ac 
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Figure 2. First example: graded Smith normal form reduction for a persistence chain complex. 



Writing the boundaries on this basis, we get a boundary module with a basis given by r\, r-i, r^, r±, r$, re, rj, 
and images in the cycle basis given by 



r\ m> tz2 — tz\ 
r 2 !->• 2 3 - tz 2 



tz^ — tz3 
T4 i — ^ tZ4 — t 2 zi 



r 5 ^ t 2 z 3 - t 3 Zl 
r 6 i-> tz 6 



o 2 

rj i — ^ t' Z5 — t Z(j 



This is the presentation map for the finitely generated lk[t]-module that represents the persistent homology 
of the sequence of spaces in Figure [l] and by computing a Smith normal form, we extract a barcode from 
the example. We shall trace the corresponding matrix and its modifications in Figure [2j In each step, we 
mark the chosen entry /y for the next clearing out. 

From this endstate, we can easily read off the barcode - especially knowing that all our operations have 
maintained the degrees of rows and columns; we have a new basis z[, z' 2 , z' 3 , z' 4 , z' 5 , z' 6 of the cycle mod- 
ule, and in this basis, the boundaries are given by tz' 2 , z' 3 ,tz' 4 ,t 3 'z' 5 ,tz' e . The resulting barcode has entries 
(1, 00), (1, 2), (2, 2), (2, 3) in dimension and entries (3, 6), (4, 5) in dimension 1. 

In above example, we can read off the map / directly since the filtration and simplicial complex is small. 
The map here is the map from the space of boundaries into the space of cycles B — > Z. By the property of 
the boundary operator, d ■ d — 0, B C Z. Hence, we can express each element of the boundary basis as a 
linear combination of elements in the cycle basis. This gives us an equation for each boundary basis element: 

hi = fij Zj 
3 

we can find fij by reducing hi with respect to Z, giving us an explicit representation of the map. 

Our second example is more abstract illustrating how the algorithm works without an explicit chain space. 

Suppose a persistence module has one presentation given by five generators x, y, z, u, v in degrees 1, 1, 2, 3, 
3, and four relations z + tx + ty, u + t 2 x + t 2 y, tv + t 2 z + t 3 y, tu + t 2 z + t 3 y. We illustrate the computation 
of its graded Smith normal form in Figure [3J From the end-state of the computation, we can read off 
a new presentation with the generators x, y' , z', u', v' and relations z' , it' ,tv,t 3 y' , where y' = y + 2x,z' = 
z + ty + tx,u' = u + t 2 y + t 2 x, v' = v — t 3 x. 

It is well worth noticing that the algorithms in existence (such as 17, 2(J 28, 3fl]) perform essentially the 
same tasks as we did using the Smith normal form to find a barcode for a homology computation; but taking 
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Figure 3. Second example: graded Smith normal form for the persistence module given 
by (x, y, z, u, v) / (z + tx + ty, u + t 2 x + t 2 y, tv + t 2 z + t y, tu + t 2 z + t 3 y). The computation 
extracts the new basis elements y' = y + 2x, z' = z + ty + tx, u' = u + t 2 y + t 2 x, v' = v~t 3 x. 
In the final matrix, the elements corresponding to finite bars are shown: two bars of length 
0, 1 of length 1 and one of length 3. The zero row corresponds to an infinite bar. 



significant shortcuts motivated by the particular cases the algorithms are built to deal with. In particular, 
most if not all algorithms only recover the pivot for the Smith normal form, and only encode the Ik [t] -module 
structure implicitly, in chosen basis orderings. 

While this algorithm is certainly not the fastest option for computing persistent homology, we believe it is 
the theoretically easiest to understand. We include a more computer science oriented presentation of the 
algorithm in Appendix [5J 



3.2. Constructions. 
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3.2.1. Presentation of Modules. Suppose M is a finitely presented module over a ring R. Then, M is given 
by some quotient R d /K for K a finitely generated submodule of R d . If R has global dimension at most 1, 
K is a free module, and thus has a finite basis expressed as elements of R d . 

Hence, to track M as a module, and to compute with elements of M, it is enough to keep track of d and of 
a basis of K. If the basis of K is maintained as a Grobner basis, then we can compute normal forms for any 
element of M that have the property that if normal forms of m, n are equal, then the elements are equal in 
the quotient module. 

Quite often, we shall meet constructions where the module M is more naturally expressed as a subquotient 
of some semantically relevant free module, possibly of higher than necessary rank. In this case, we shall 
track two sets of data to enable computation in M; a basis of the relations submodule K, and an extension 
of that basis to a basis of the generators submodule G, resulting in two bases that track the behaviour of 
M = G/K. Again, as long as K has a basis maintained as a Grobner basis, it is enough to compute this 
normal form to compare elements. 

In particular, this is the case where the module M is the homology of some chain complex. There, the free 
module is the module of chains, and the submodules G and K correspond to the modules of cycles and 
boundaries respectively. 

We observe that this reflects practice among algorithms for persistent homology: to compute the persistent 
homology with a basis of a filtered simplicial complex, we maintain a cycle basis and a boundary basis, 
reducing modulo the boundary at each step to find out if the boundary da of a newly introduced simplex a 
is already a boundary or not. 

3.2.2. Useful Forms. We have already seen the introduction of a Smith normal form to the graded k[t]- 
module context in Section [3.1.31 The Smith normal form transforms a matrix representing a map between 
free modules by basis changes in the source and target modules until the matrix is diagonal; all while 
respecting the grading of the modules. 

Another fundamentally useful form that we will be using a lot is the Row or Column echelon form of 
a matrix. For a matrix F representing a map / : M — » N between free graded k[i]-modules, a row echelon 
form changes basis for N to eliminate all redundancies in the information contained in the matrix, while a 
column echelon form changes basis for M to achieve the same goal. Certainly, a Smith normal form is an 
echelon form, but the converse does not necessarily hold true. 

The properties most interesting to us for an echelon form - and we shall state these for the row echelon form; 
column echelon form holds mutatis mutandis — are: 

• All completely zero rows are at the bottom of the matrix 

• The leading coefficient of a non-zero row (the pivot of that row) is strictly to the right of the leading 
coefficient of the row above it 

• All entries in a column below a leading coefficient of some row are all zero. 

The power of the echelon forms come in what they imply when you use a matrix to pick out a basis of a 
submodule: if the basis elements of a submodule of a free module are in the rows of a matrix, and this matrix 
is placed on row echelon form, then several problems concerning module membership and coefficient choice 
become easy to work on. 

To determine if an element is in the submodule, add the element to the bottom row of the matrix, and put 
the matrix on row echelon form again. This, by force, will eliminate any entries in pivot columns of the new 
row; and modify elements to the right of the pivots that can reduce. Eliminating all column entries that 
occur as pivots of the previous basis puts the new element on a normal form with respect to the basis - 
and this normal form is sufficient to determine equality modulo the submodule represented by the basis. 

3.3. Representation of Homomorphisms. Cohen-Steiner, Edelsbrunner, Harer, and Morozov [IH de- 
scribe a method for computing the kernel, image, and cokernel persistence for the case of what they call 
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compatible nitrations. For two filtrations, / over a space X and g over a space Y, the two nitrations are 
compatible if X C Y and the restriction of g to X is equal to /. 

The algebraic equivalent (and generalization) of this is the notion of graded maps. 
3.3.1. Dictionary of Operations. 

3.4. Two Modules. When given pairs of persistence modules, we must be able to represent the maps 
between them, assumin g of course, they follow the assumptions in Section 13.31 This material follows the 
treatment by Eisenbud 2l|, Chapter 15], specialized to the case of persistence modules. We begin with 
presentations for modules P and Q: 







^G P - 

: f\a t 



F P 
-> Fo - 



-> P 



l-'Q 



-> 



-> 



Such a presentation exhibits P and Q as quotients of free modules, i.e. P = Fp/ipGp and Q = Fq/zqCq. 
We call Fp and Fq modules of generators of P and Q respectively, and Gp and Gq modules of relations. 

Over nice enough rings^, both Gp and Gq are free for all modules P. These rings include all Euclidean 
domains. 

Now, an arbitrary map P — > Q can be represented by a map from the generators of P to the generators of 
Q, or in other words a map Fp Fq. For such a map between the generator modules to represent a map 
between the quotient modules it needs to obey one condition: ip ■ ipGp C %q ■ Gq. In other words, any vector 
in Fp that represents a relation in P, and therefore represents the 0-element in P, has to be mapped to a 
relation in Q so that its image still represent the 0-element. Subject to this condition, any map between the 
free modules Fp and Fq will represent a map between P and Q. 

Among the core messages of this paper is the recognition that by maintaining Fp and Gp we can compute 
just about anything to do with persistence modules. In particular, if Fp is kept in a nice form - so that all 
the elements of Gp form a tagged sub-basis of the basis of Pp, and all basis elements of Pp not already in 
the Gp basis are kept reduced with respect to Gp, then questions about normal forms for elements of P, 
equality in P, and many other relevant questions become easy to handle algorithmically. 

Throughout the expositions below, it is universally helpful to keep any module of relations on a row reduced 
echelon form and any module of generators reduced and augmented by the relations. By doing this, we can 
consider the effects on single basis elements of the various constructions, and thereby the effects on a barcode 
presentation of a module. 

This kind of nested module presentation is to some extent inherent in the persistence barcode: for the 
barcode to be accessible, the bases for the persistence module expressed are kept in a format where any 
finite barcode gives rise to a pair of basis elements - one of Gp that maps onto a basis element in Pp, 
with some grading shift corresponding to the length of the barcode. By keeping the bases for Fp and Gp 
synchronized so that all basis elements of Gp form halves of such pairs, the barcode becomes immediately 
connected to the basis presentation. 

To accomplish any useful computation with these representations of persistence modules, we need to be able 
to have some basic tools or manipulate them. In all cases, this will be done in terms of some submod- 
ule of a joint quotient module. That is, by considering both modules simultaneously and extracting the 
relevant information we will end up in a situation where we input presentations and as output get presenta- 
tions, thus allowing these tools to be composed together algorithmically. In addition, we relate the matrix 
representations more familiar to persistence literature to each case. 



Global dimension less than or equal to 1, which includes all PIDs and therefore all Euclidean domains 
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3.4.1. Chains or coefficients. It is worth pointing out that depending on the application domain, two different 
representations for the presentation may both be worth while. One way to represent the generators and 
relations, especially in an algebraic topological setting, is by giving explicit chains for both in a chain module 
containing both. This way, the modules Gp and Fp are represented, essentially, by matrices representing 
the maps Gp — > C* and Fp — > C* for the chain module C*. 

In this setting, algebraic manipulation on generators and relations can be performed separately from each 
other. In particular, the relations module can be set on a row echelon form independently of the generators 
module. This is also the format emitted by existing algorithms and implementations. One major drawback 
of this representation is that the map Gp — > C* has to be pulled back to a map Gp — > Fp before a barcode 
can be computed. 

Alternatively, the two maps Gp — > Fp and Fp — > C* can be maintained separately. In this setting, the 
components of a barcode for P are easy to extract; and it seems likely that the size of the matrix Gp — > Fp 
will be significantly smaller than a matrix Gp — > C*. However, each basis change in Fp comes with the 
requirement to modify two different matrices - for each operation on the matrix Fp — > C», the inverse 
operation has to be performed, transposed, on the matrix for Gp — > Fp. 



3.4.2. Barcodes as module presentations. For the particular case of a persistence barcode, the setting at 
hand is very specific. The modules are N-graded modules over k[t], and in the presentation of a persistence 
module P as 

O^Gp^Fp^P^O 

we maintain Gp and Fp as free submodules of a global chain module, with a global specific basis given by 
the simplices in the underlying simplicial complex. In this setting, we further maintain a few invariants, in 
particular we ensure that Gp and Fp are the results of the appropriate basis changes S and T from the 
Smith normal form in order to guarantee that ip is the corresponding diagonal matrix representation of the 
Smith normal form. 

Hence, from an arbitrary presentation of a finitely presented persistence module 

O^Gp^Fp^P^O 

we get the barcode presentation by writing tp = SipT and then replacing this presentation above with 

-> T^Gp ^ FpS- 1 -^P^O 

One concrete benefit of this approach is that in the new bases for T _1 Gp and FpS , we can match up 
basis elements between the relations and the generators into birth/death pairs of barcode intervals, and the 
corresponding diagonal entries in tp are exactly on the form kt a yielding a bar of length a in the barcode. 



3.5. Direct Sum. The first construction is the direct sum of P and Q. This is the simplest and most basic 
construction, but will be used in all the other constructions. We construct it by taking the direct sum across 
the relations and free part. 

(ip 0\ (pp 0\ 

^ G P © G Q -i Fp © F Q -i P © Q -> 

This is still a short exact sequence since with a direct sum the maps can be composed as a direct sum of the 
component maps. 

In the case where the persistence module is represented as a matrix encoding a cycle basis and another 
matrix encoding a boundary operator, the new representation is simply a block matrix with the component 
matrices along the diagonal. Consider the matrix equations 



Rp = DpCp 
Rq = DqCq 
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Rp 


" 




'Dp " 




C P " 





Rq_ 




Dq 




. C Q. 



The direct sum is given by 
(1) 



Assuming that the two inclusion maps i p and iq are already on Smith normal form, the Smith normal form 
of the presentation map ip ®iq is given by the block matrix 

i P 

%Q 



3.6. Image. Given a map, we will want to construct the image of the map between the two modules. This 
is a relatively straightforward construction. We first compute the image of the map between the free parts: 
img</3 : img(i 7 p — > Fq). The finite presentation is then 

-> G Q G Q + img^ -> img(P 4g)^0 

Note that the resulting image generators should be reduced with respect to the relations to ensure that they 
are representative generating elements. 

It is sufficient to only consider Gq since we require tp ■ ipGp C iq ■ Gq. If tp is an inclusion, this implies 
that Gp C Gq. This is a general statement which implies that if we have an inclusion, generators will map 
to generators and boundaries map to boundaries. 



In the context of this in topological terms as presented in [la], we choose a basis for Fp and Gp. By inclusion, 
img<^ = Fp. Since Gp C Gq, Gp is presented in a compatible basis for F p and so we simply extend this 
basis to Gq. Algorithmically this means finding computing the persistent homology of P then adding the 
simplices which are in Q and not in P if they add any relations (positive simplices are ignored as these would 
be in Fq). 

If we are starting with both module presentation on a barcode form, with the inclusion maps ip and iq both 
on Smith normal form, then we can approach a Smith normal form for the presentation of the image by first 
computing the sequence of all <p(fj) for fj the basis of Fp used. 



Now consider the block matrix 



/ •¥>(/!).. 
• ■■<p(fr) ■ ■ 

■ ■■■ip(gi) ■ 
V ■■■ip(g s )- 



\ 



Ir 



I s 



J 



with all of the images expressed in some basis of Fq, and the right hand blocks consisting of an identity 



matrix in order to track the changes that will happen, 
we will find a matrix on the shape 



* 
.4 



Putting this matrix of a row reduced echelon form 



where the blocks (AB) carries the coefficients expressing the elements ip(gj) in terms of the collection of 
elements ip{fj) and ip(gj). In other words, this part is a presentation matrix for the map Gq — > Gq +img ip 
that we need. Putting this on Smith normal form yields barcodes for the image module. 

3.7. Cokernel. To compute a cokernel, we merely include the images of the generators in Fp among the 
relations for the cokernel as a quotient module of Fq. 



O^Gc 



F P F Q -> cokcr(P 4(5)^0 
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To see why this is correct, first consider the cokernel of the free modules. In this presentation if it is in the 
image, since we have listed it among the relations, it maps to in the quotient module. One imaginable 
complication would be if Fp has a basis element representing in P, but mapping to a basis element in 
Fq which is not killed by any element in Gq. For a constellation like this, we would in the presentation 
above erroneously kill the generator in Fq by the presence in Fp; but the constellation is impossible, since 
(p ■ i P ■ G P C i Q ■ Gq. 

In terms of algorithms [15j . we again compute a basis for Fp and Gp just as in the image case. We now 
only consider cycles and relations which do not appear in P. 

3.8. Kernel. In order to construct the kernel for arbitrary finitely presented modules, we first describe how 
to construct it for free modules. Suppose we have a map between free modules ip : F — > G. Then the kernel 
of this map is the free submodule of F consisting of elements that map to in G. This can be computed 
by the usual Gaussian elimination approach: we can write down a matrix with rows in G © F, each row a 
pair (ip(f), /) for / a basis element of F. After putting this matrix on a row-reduced echelon form, some of 
the rows will have entries only in the F part, all entries in the G summand having reduced to 0. Projecting 
these rows back to F gives us a basis for the kernel. 

The kernel of a map between finitely presented modules has to be constructed in a two-step process. A 
generator of the kernel of P — > Q is an element of Fp such that its image in Fq is also the image of 
something in Gq . This is to say that the free module of generators of ker(P ^> Q) is given by the kernel of 
the map between free modules Fp © Gq — > Fq given by (f,g) H> (f(f) — *Q (<?))• 

This takes care of the module Fk of generators of ker(P Q). For a complete presentation we also need 
the relations of this module. Since the kernel module is a submodule of P, the relations are the restriction 
of the relations in P to the kernel K. These can be computed by recognizing that they are given precisely 
by the elements in Fk such that their image of Fk in Fp coincides with images of elements in Gp. Thus, 
we setup another map Fk © Gp —¥ Fp, defined by (/, g) i-> fp — ip(g), and compute the kernel of this map. 
This kernel is the module of relations we need. 

For this case, we cannot simply choose a basis and extend it in the algorithmics (as the more complicated 
structure described above shows). There is no immediate way to take short cuts around the Smith normal 
form computation either. 

Notice that if elements of the kernel are computed in some alternative way, then the following lemma may 
be useful: 

Lemma 8. If M is a free k[t] -module and cf> : N — > M is k[t] -linear, then if tz € ker^, then z € ker</>. 

Proof. Since M is free, the only way that tw = in M is if w = 0. Since tz £ ker (j>, we know that d(tz) = 0. 
But by k[i]-linearity, d(tz) = td(z). Hence, td(z) = 0, so d(z) = follows. □ 

In particular, if (C, d : C — > C) is a free k[i]-linear chain complex, then if tz is a cycle, then so is z. 

3.9. Free Pullback. Before continuing into new constructions, we recount a useful construction on free 
modules: the free pullback. Given two maps, / and g in the following diagram: 

.4 

/ 



we construct the pullback P such that the following diagram commutes 
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P 



> A 



D 



-> C 



To compute P, we must set up a kernel computation: namely find a basis of ker(/ < 
construct the following matrix 



-g). Graphically, we 




We perform a column reduction on the matrix [F\ — G] (on the right) to compute a basis for the kernel, 
keep track of the basis we perform the same operations on the matrix on the right. 



To 



HH1 



Ml 



Once computed, the matrix on the right has some zero columns, shown above in blue. The corresponding 
columns in the matrix on the right, represent a basis for the kernel. The projection maps are given by taking 
the appropriate part of the basis: the top part (the darker blue) represents the projection onto the domain 
of F and the lighter blue in the lower part represents the projection onto the domain of G. Note that the 
elements in the top right quadrant, represent the kernel of F and the elements with only non-zero entries in 
the lower half represent a basis for the kernel of G. 



With this construction we can restate the kernel construction in Section 
the kernel Fk are given by the pullback 

Fr ' ' > Fp 



graphically. The generators of 



G c 



^Fr 



The relations Gk are given by a second pullback 



G 



K 



> Gi 



Fk — ^— > Fp 



giving the presentation Gk — ^> Fk — > ker (<£>). 
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3.10. Pullback. Given module presentations 



- 


■» G P 




^P^O 


- 


^G Q 


^F Q 




- 


-> Gr 


^F R 


^R^O 


and maps P A- R and Q A R given as 


maps 


F P A 


F R and Ft 



ib 

— > Fp such that (p(Gp) C Gp and 
^{Gq) C Gp, we can dehne the pullback. This is the submodule of P © Q such that for any element (p, q), 

ipp = ipq. 

Pullbacks are useful when dealing with modules - they can express very many algebraic concepts neatly and 
concisely, such as inverse images of maps, fibre products, kernels, and so on. 

To compute the pullback, we set this up analogously as with the free pullback, as a kernel computation. The 
pullback is a submodule of P © Q, and is carved out exactly as the kernel of the map P © Q —> R given by 
(p, q) h-> ipp — tpq. Since this is the kernel of a module map, we can compute this using the construction of 
Section [ 







-> Gi 



G c 



<p\g p — '4>\g c 



-> Gi 



-> F P ®F Q - 
(p—ip 
^F R 



-> P( 



-> 



f-g 



-> J? 



-> 



Expanding the construction, we compute the following two free pullbacks 



Fpb F > Fp 



F c 



^Fp 



Grb 

"V" 

Fpb 



> G ? 



-> Fp © Fq 



Here we begin to see that with presentation algorithms can be broken down into a few key constructions. 



3.11. Pushout. Pushouts are the dual construction to pullbacks. Some of their most important uses is to 
glue together things that overlap slightly - producing almost but not quite a direct sum. 



Given module presentations 



O^Gp 
O^Gq 
O^Gr 



IP, 


F P 


Pp 


p - 


+ 


ia 


Fq 


po v 


Q- 


» 


in ^ 


Fr 


PR 


R - 


* 



and maps R — > P and F — > Q given as maps Fp — ► Fp and Fp — !> Fq such that y(Gp) C Gp and 
iP(Gr) C Gq, we can define the pushout. 

This is the cokernel of the map F — >• F x Q given by r i— ► (yr, —tpr). This way, any image of an element F 
in either P or Q can "move across the x" in the pushout module P XrQ. 



3.12. Tensor products. The tensor product, and its various associated constructions, are easiest to de- 
scribe if we focus on using the presentation map G — > F instead of explicit chain representations for the 
relations. We assume that readers have seen tensor products at least in the context of vector spaces and we 
recommend Eisenbud [2l| as a reference. 



We fix a basis Bp for Fp and Bq for Fq. A basis for M © N is given by Bp © Bq. 
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Relations are all generated by elements on the form r ® g or g ®r where r is a relation and g is a generator. 
Theorem 9. Assume that the underlying coefficient ring is a graded PID. 

Suppose P is presented by i : Gp — > Fp given on Smith normal form, with Bp the corresponding basis of 
Fp, and Q is presented by j : Gq — > Fq also given on Smith normal form, with Bq its corresponding basis 
°f Fq. 

Then a Smith normal form presentation of P ® Q is given by a basis Bp x Bq for the generating module, 
and a basis element in the relations module for each p ® q for p G Bp, q G Bq where at least one of ap and 
(3q is a relation. Let 7 be the generator of the ideal (a,/?). Then jp®q is the sole relation influencing p ® q 
in the Smith normal form presentation of P (£> Q. 

Proof. First, if all relations for both P and Q are given as Smith normal forms, any relation on the form 
ap £g> q for p, q images of the relations basis are going to be products of basis elements; and thus not expand 
bi-linearly. Hence, each relation induced by the relations construction above is already a coefficient times a 
basis element. 

To ascertain a graded Smith normal form, we also need to verify that each basis element occurs in only one 
relation. This is not immediately guaranteed from the construction above - there is one small correction 
step needed. Certainly, if ap G iGp and j3q G ]Gq, then jp ® q has to be a relation. However, if both 
and Pq 7^ 0, then both ap C5> q and fip ® q will show up from the construction above; for these cases, 
replacing the two candidates by their common generator jp ® q is needed to construct a Smith normal form 
presentation. □ 

We may note that the tensor product of the presentation is longer than the presentation suggested above: 
it would be 

-)• G P <g> Gq ->• G P ® Fq © Fp <g> Gq ->- Fp ® Fq -)• P ® Q ->• 

This is a free resolution oi P ® Q; but not a minimal one — the method in the proof of Theorem [9] tells 
us exactly how to reduce away the redundancy represented by the syzygies in Gp ® Gq to get a minimal 
presentation — of length corresponding to the homological dimension of the module category. 

For an arbitrary size tensor product P\ ® ■ ■ ■ ® P n with Pj having relations Gj , generators Fj , Smith normal 
form presentation map ij, and a basis of Fj denoted by Bj, a basis for the generators is given by n-tuples 
in Bi X • • • X B n , and relations are given from the generator of the ideal generated by all the coefficients of 
relations of factor basis elements of the relation. 

We can summarize the algorithm above as follows: 

Tensor Product Algorithm 

Input: Two module presentations P and Q with reduced bases in Smith Normal Form 

(1) The generators are the tensor product of the generators: F p ®q = F P ® Fq — generating elements 
are of the form (fp(i), /qO)) for all generator elements in P and Q 

(2) Create a non-minimal list of relations by pairing relations and generators — Gp^Q = Fp ® Gq © 
F Q ®G P 

(3) Create a minimal representation of relations — reduce the relations modulo Gp ® Gq — consider 
the pairs (<7 p ,/q) and (/ p ,<?q) such that g p — > f p and gQ — > /q and only keep the relation which 
occurs first. 

Output: A presentation of the tensor product 

We give an example of this process in Section [33221 

3.12.1. Interpretations of tensor products. The tensor product described above is the tensor product M ®^t t j 
N, in which the t parameter can be freely moved "across" the <g) sign. One can consider this to be a tensor- 
product of persistence modules, in which time steps for the conglomerate by stepping time arbitrarily on 
one side or the other. 
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Algorithm 1 Tensor Product Algorithm 

1: Input: Two module presentations P and Q with reduced bases in SNF: 
Gp^Fp^P 
Gq^Fq^Q 

2: Compute generators: F P ® Q = F P ® F Q = (/p(i),/q(j))V/ p (j) G F p , /q(j) G Fq 

3: Compute non-minimal relations: Gp®q — Fp ® Gq © Fq ® Gp = (fp(i), 9q{j)) © 

(ffp(*), /qW) V/ P (i) G F P , / Q (£) g F q , 3p (j) g Gp, / Q (fc) G F Q , 
4: For pairs (g p , /q) and (/ p ,5q) such that g p — > f p and <jq — ► /q, keep the relation which occurs first. 
5: Output: Relations, generators and a map: Gp^Q — » Fp^Q 



There are more tensor products we can construct: by forgetting about the persistence structure on one or the 
other side - or by having a persistence module M and a Ik-vector space N - we can produce tensor products 
M ®^ AT. For these, the t action does not move fluidly between different factors of the tensor product; but 
we can still recover a lk[i]-module structure on the tensor product. 

In fact, if we have persistence modules M,N, then M ®^ N has two different k[£]-module structures; these 
are defined on the basis elements by 

t(m ® n) = tm © n t(m ® n) — m ®tn 

and extended by linearity to the entire tensor product. We'll explain the first of these structures - the 
second one is completely analogous. In this case, different persistence parameter values for N give completely 
different components of the product: stepping with time only moves up on the M side of the factor. 

These constructions correspond to the different dualizing functors introduced by de Silva, Morozov, and 
Vejdemo- Johansson [l7| . 

3.12.2. Relations to horn modules. Suppose M is a persistence module. Then M* — hom(M, k[f]) is a 
persistence module too, whose degree d elements are maps M — > k[t] that change degree by d; so that if 
1 771 1 = n and |/| = d then |/(tt7,)| = n + d. 

This module has the structure of a graded k[i]-module by the rule (tf)(m) = f(tm), and we can easily verify 
that this (tf) has an increased degree as defined above. 

Given a basis B = {mi, . . . , rrik} of M, there is a dual basis B* = {m\, . . . , of M* where m*(m,j) is the 
Dirac delta function 8ij . 

We consider the module hom(M, N) for persistence modules M, N consisting of maps / : M — > N of 
arbitrary degree, just as defined above: / has degree \d\ if when |m| = n then |/(m)| = d + n. 

Example: Suppose M is given by two generators x, y in degrees 1 and 2, and the relations t 3 x, t 4 y. Suppose 
that N is given by three generators u, v, w in degrees 1,1,2, and with relations t 2 u, tv, t A w. The dual module 
M* has generators x*,y* : M -> k[t] defined by x*(x) = l,x*(y) = 0,y*(x) = 0,y*(y) = 1. Thus, the 
generators have degrees —1 and —2 respectively. We can establish by inspection that t 3 x* — and t 4 y* = 0, 
producing the relations of M*. 

A generating set for M* © N is given by all pairs of elements: x* ® u, x* <E> v, x* ® w, y* ® u, y* ®v,y*® w. 
Relations are all pairs of generator with relation, thus given by the following list: 

t 3 x* <g> u t 3 x* ® v t 3 x* <g) w t 4 y* ® u t 4 y* © v t 4 y* ® w 

x* (8) < 2 m x* ®tv x* ® t 4 w y* ® t 2 u y* ®tv y* ® t 4 w 

Out of each column, the lowest degree of t remains as the actual relation, giving us the final list t 2 ■ x* ® 

u,t ■ x* ® v,t 3 ■ x* ®w,t 2 ■ y* ® u,t ■ y* ® v, t 4 y* ® w. 
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One persistence homomorphism / : M — > N is given by f(x) = u + v, f(y) = tu + w. In the presentation 
above, this homomorphism is given by the linear combination 

f = x*®u + x*®v + t-y*®u + y*®w 

and by computing degrees of each term, summing degrees for the basis elements in each pair, we can see 
that / has degree as a homomorphism. 

3.13. Exterior powers. The exterior power A 2 (M) is given by M (g)M/ ~ where a<S)b ~ —6(8) a. We get it 
by adding / ® g + g (g> / for all generators /, g of M to the relations module of the presentation of the tensor 
product M ® M. 

We can choose a different representation and handle the exterior powers. The appropriate representation 
here is to assign as a generator of the exterior power A m (M) a set of m distinct elements from the component 
basis elements, paired with a permutation a e S m tracking the "original" order of the basis elements. This 
permutation is important for sign handling when comparing different elements - (M, a) = (M, r) iff err -1 is 
an even permutation; if ot -1 is an odd permutation, then (M, a) — — (M, r). 

3.13.1. Discrimination with exterior powers. The exterior power construction has some discriminating power 
between possible barcodes, even without computing a Smith normal form directly; just computing the exterior 
power can in some cases distinguish between different possibilities for barcodes. 

As an example, consider a persistence module M with two finite barcodes. It has generators in degrees 
1, 2 and relations in degrees 5, 10. Just from degree reasons, we know that the possible barcodes are either 
(1,5); (2, 10) or (1, 10); (2, 5). 

Call the generators x, y; the exterior power A 2 M has as its generating set the set of all cardinality 2 subsets 
of {x, y}, which is to say the only basis element is x A y. 

As for the relations, a relation at d x + bt e y gives rise to relations at d (x A y) and bt e (x A y); which together 
have as their basis ct^ (x A y) where ct^ is the higher-degree choice out of at d and bt e . The higher degree 
choice is really just the gcd and is due to the nature of basis changes in these modules; if x exists earlier 
than y, then y can be replaced by y — Xt k x, but x cannot be replaced by x — \t~ k y. Since the relation is 
homogeneous, the earlier existing basis element is the one with the higher degree power of t associated. 

The presentation matrix for M will be 

fat 4 ct 9 \ 
[bt 3 dt s ) 

This yields a presentation matrix for A 2 M that has entries in degree order 

(-at 4 \ 
at 4 
bt 3 
-bt 3 
-ct 9 
ct 9 
dt 8 
\-dt*J 

with Smith normal form dependent on which of a, b, c, d are non-zero. In particular, if b ^ 0, then t 3 x A y is 
a relation, and if b = then t 4 x A y is a relation. 

Consider now the two cases we gave here. In the first option, our relations are t x, t s y, and in the second 
they are t 9 x, t 3 y. If the relation in the exterior square is t 4 x A y or t 3 x A y, this tells us which of the two 
options applies to our case. 

This agrees with what we would have expected from a computation of a barcode for M before we computed 
the exterior algebra. Indeed, if b ^ 0, then the Smith normal form computation of the presentation for M 
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would have started by a change of basis, replacing y by y' = y + a/btx. As a result, the relation transforms 
into bt 3 y', giving a barcode (2, 5) which settles the entire module structure. 

3.14. Symmetric powers. The symmetric power S 2 (M) is given by M® Mj ~ where a®b ~ b ® a. We 
get it by adding f ® g — g ® f for all generators f,g of M to the relations module of the presentation of the 
tensor product M (g> M. 

Alternatively, one may represent the symmetric power S m {M) with a basis corresponding to multisets of 
weight m of basis elements. The same construction holds for finding relations as for the tensor products - 
with the generating element of the coefficients that kill each factor basis element in the multiset killing the 
compound basis element. 

As opposed to exterior powers and the example in Section 13.13.11 symmetric powers force commutativity 
instead of anti-commutativity. Where A k (M) eventually vanishes for finite dimensional M, S k (M) is finite 
dimensional and non-trivial whenever M is. 



4. Applications 



Although our primary goal was to describe a general framework for manipulating persistence modules in a 
common framework, there are several application areas that have motivated the constructions in this paper 
and show great promise in applicability for the algorithms and ideas present. We shall discuss some of these 
here. 



4.1. Torsion complexes and persistent relative homology. Consider a persistent chain module with 
torsion - a [k[i]-module with some non-trivial relations, and with an endomorphism d of degree that squares 
to 0. Such modules arise in various contexts - as part of stepwise approximation in spectral sequence 
approaches to persistent homology, as well as when computing relative homology in a persistent setting. In 
the former, torsion comes from the fact that we are dealing with persistence modules and in the latter comes 
from the quotienting at a space level. 

Now, the persistent homology of this directed system C* is simply homology with coefficients in k[t] of 
the torsion persistence module. In other words, we can represent torsion at the chain level and compute 
homology as we normally would using the techniques described earlier in this paper; we have all the tools 
necessary for the computation of persistent homology of a torsion chain complex. For a complete example 
of computing persistent relative homology see Appendix [B] 



4.2. Spectral sequences. This paper was originally motivated by our work on spectral sequences in per- 
sistent homology [26(. For our interest in spectral sequences where each entry is a persistence module, all 
fundamental algebraic constructions need to be accessible - to compute subsequent pages of the spectral se- 
quence, one repeatedly computes homologies locally at each point of a boundary map that successively 
approximates an end-result closer and closer. 

For the first step, going from E° to E 1 , the modules can usually be assumed to be free modules, and the 
justifications by Zomorodian and Carlsson (35j are enough to assert computability; but even for spectral 
sequences going from E 1 to E 2 and beyond, a firm grasp of the computation of kernels and cokernels of 
persistence modules with possible torsion is needed. 

With the algorithms stated here, we are able to produce code for computing spectral sequences with persis- 
tence modules. These spectral sequences could be bases for several interesting directions for further study: 
parallelization by using a spectral sequence to merge local information into global, as described in [26| . or 
computing homology of the total space of a combinatorial fibration using fibres and base space by adapting 



the Serre spectral sequence as described by Goerss and Jardine [24J are two approaches that come to mind 



4.3. Unordered persistence computation. Many of the most efficient algorithms for constructing filtered 
complexes in the first place, such as the algorithms proposed by Zomorodian and Carlsson [35j, produce sim- 
plices in some order compatible with the inclusion ordering on the simplices, but not necessarily compatible 
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with the chosen filtration ordering. As an alternative to re-sorting the simplices this algebraic approach 
offers us algorithms to compute with an essentially unordered simplex sequence, performing reductions as 
needed, based on the new simplices acquired. 

If the critical values of the filtration function are not known before hand, everything in this section can 
still be performed as long as one is prepared to adjust degrees for all elements in play along the way; for 
exposition, we shall assume that each simplex has a global N-graded degree associated to it. 

The classical persistence algorithm as described by Edelsbrunner, Letscher, and Zomorodian [20| and by 
Zomorodian and Carlsson [35j works by computing a Gaussian elimination on the boundary matrix, respect- 
ing the grading on the rows and columns involved by the choice of ordering. 

When introducing a new boundary basis element, this will affect the current state of the Gaussian elimination 
running under the hood by first reducing the new basis element with everything that was already there - 
by filtration value - when the new element showed up, and then using the thus reduced element to further 
reduce any later additions to the boundary basis. This effect of the grading is transparent enough that we 
can fit it in a variation of the original persistence algorithm. 

Adding a new simplex a, we compute da, and reduce da by all the boundaries that were added at earlier 
filtration values than that of a. If the result is 0, then clearly, da was a boundary, and so a is the leading 
simplex of a cycle. 

If da does not vanish when reduced modulo the current boundary basis, then the leading simplex of the 
reduction of da is some simplex r associated to the cycle basis. Thus, a is a pre-boundary, bounding the 
cycle that has r as its leading term. We introduce the interval (/(r), f{a)) as a new interval. 

Invariant of the persistence algorithm: at any stage, the simplices are paired up in a partial matching such 
that if (r, a) is a pair, then da reduces modulo all the boundaries existing at the time of the introduction of 
a to a chain with leading simplex r. 

This invariant can change with the introduction of some new simplex v in the following ways: 

• dv reduces modulo earlier boundaries to 0. Thus, v starts a new cycle. 

• dv reduces modulo earlier boundaries to some cycle with leading simplex 9. 

— If 6 is the leading simplex of some later boundary basis element, linked to a generating pre- 
boundary simplex r, then there is some pair (6, r). If f(v) < /(t), then this pair is modified to 
{9,v) and the boundary linked to r is further reduced using dv. 

This may well lead to cascading changes, as the resulting new boundary chains in turn reduce 
later and later chains, and replace more and more paired simplices. 

— Otherwise, 9 is an unpaired cycle simplex (since the ordering is compatible with inclusion of 
simplices), and therefore creates a new simplex pairing and a new interval in the barcode. 

4.3.1. Example. Take the 2-simplex spanned by 1,2,3 and introduce the simplices in the following order, 
with the filtration values given as well: 



Simplex 


1 2 


12 


3 


13 


23 


123 


Filtration value 


1 4 


6 


2 


3 


5 


7 



At the start of the algorithm, nothing is out of the ordinary - as long as filtration values increase monoton- 
ically, this works exactly as the classical persistence algorithm would. Hence, after consuming 3 simplices, 
we have a cycles basis consisting of 1, and a boundary basis consisting of 2 — 1, and a pairing of (2, 12). 

Then arrives the simplex 3. We compute 9(3) = 0, and add 3 to the cycles basis. 

Next comes the simplex 13. The boundary is 3 — 1, and we pair it with 3, yielding the pairings (2, 12), (3, 13). 
At this stage, we have a cycles basis of 1 and a boundary basis of 2 — 1, 3 — 1. It is worth noting at this 
point that everything involved knows of its filtration value, and therefore of its own grading. 

Now, the simplex 23 arrives. The 1-simplices are ordered by their filtration values as /(13) < /(23) < /(12), 
which concretely means that 9(23) = 3 — 2 only gets reduced by 3 — 1, to the value — 2 + 1 = —9(12). 
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Further reductions are impossible with the older parts of the boundary basis, and thus this simplex should 
have originally paired with its leading element. We adjust the pairing (2,12) to (2,23). Furthermore, the 
previously not entirely reduced element of the boundary basis given by 9(12) = 2 — 1 now is reducible with 
this new simplex information; at which point 12 transforms from a pre-boundary to a cycle. 

Finally, the introduction of 123 kills the cycle 12 as expected. 



5. Conclusions & Future Work 



The primary motivation is to provide a general framework for more complex operations with persistence 
modules — hopefully reducing the need to design new algorithms from scratch for new constructions. The 
trade-off is that these are not optimal in terms of running time or space but provide the flexibility for 
experimentation. 

We have demonstrated how techniques from commutative and computational algebra influence the study and 
the computation of persistent homology. We have shown how the classical concept of Smith normal form can 
be adapted to graded modules over a graded PID, and how this adaptation leads to concrete algorithms for 
computing with persistent barcodes, and maps between persistence modules. The common representation we 
use is the presentation of a finitely generated module. By making the representation explicit, we are able to 
state algorithms which take as input presentations and output presentations, allowing for easy composition 
of different constructions. 

It would be interesting to extend the understanding of the interplay between algebra and topology - un- 
derstand how different categories of modules, different representation theories affect the computation of 
persistence. We have already seen three categories that are useful for different approaches to persistence 
- modules over k[t], modules over kQ for a quiver of type A n . Furthermore, the representation theory of 



k[s, t] and the classification by Gabriel 22J of representations into tame and wild types were instrumental to 



proving the limitations of multidimensional persistence. 
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Appendix A. Algorithm for Computing Presentations of Persistence Modules 

In this appendix we give a more explicit algorithmic description of Section 13.1.31 We begin by describing 
computing standard persistent homology. 

Given a filtered simplicial complex, we write down its graded boundary operator. The added needed data 
structure which is a data structure which returns the relative grading between two simplices. If we have a 
discrete parameter filtration (or if we discretize a continuous parameter), then the relative grading is the 
time difference between when the simplices appear. 

Therefore we can write down the boundary operator matrix using the standard formula 



where the ?Vs are the vertices of the simplex a. 

For each non-zero entry, we multiply it by the relative grading of the two corresponding simplices. Since we 
are in chain space, each row/column pair corresponds to two simplices. We adopt the following conventions: 
we sort the row space such that each column of the matrix has decreasing powers of t. This is equivalent to 
sorting the simplices in order of appearance in the filtration. We sort the columns, such that each row has 
increasing powers of t. This also corresponds to sorting the columns in terms order of appearance. 

We now introduce the standard reduction algorithm. We use the following notation 

• F - a graded m x n matrix such that columns have decreasing powers of t and rows have increasing 
powers of t 

• Fi - i-th column vector 

• Fij - (i,j)-th entry in F 

• pivot s(-) - a lookup table where we store pivot elements, it returns the column index which contains 
the pivot element of that row; if empty returns 

• diff(-, •) - a function which return the relative grading difference between the two inputs 

• C - matrix which stores chains 

• Z - kernel basis of F (cycle basis) 

• B - image basis of F (boundary basis) 

When computing persistence we set F to the graded boundary matrix d. 

Note that this gives the basis elements of Z and B in order of appearance. Each element is represented as a 
graded chain. As in the example in Section [3. 1.3[ we now must write down the map B — > Z. To obtain the 
barcode, we must reduce the map to SNF. 

The reduction algorithm of the map is almost identical to the algorithm above except that once a pivot is 
found, we reduce the rows above it and we do not need to store the change of basis if we only care about 
the barcode. We use the same notation as above, except F : B — > Z and since we will need row operations 
as well, F l represents the z-th row of F. 

To put this into context of the algorithm described in As above it computes cycles in order of appearance and 
stores them by only keeping track of the element with the lowest power of t. When computing boundaries, 
it only considers the elements corresponding to elements in the cycle basis. Therefore, the pivot already 
corresponds to an entry in the SNF, since it represents a map from a boundary element to a cycle (represented 
by one element — the last simplex which "created" the cycle). While this is more efficient, it makes it more 
difficult to generalize since we do not always have the freedom to choose such nice bases . 

To generalize the above algorithm, we must have the ability to find relative gradings between chains. In the 
above case, this can be done by mapping back to simplices. In general with more complicated operations, 



(2) 




2<> 
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Algorithm 2 Reduction Algorithm 

1: Input: F = d 

2: Initialize pivots = 0, C = I (n x n identity matrix ),Z = 0, and B = 

3: for i=l,. . . ,n do 

4: while Ft ^ do 

5: Set j «— lowest non-zero entry (smallest power of t) 

6: if pivots(j) 7^ then 

7: X = dif f {pivot s{j),j) 

8: Ft F £ I ■^i > pivots(j)^)-^pivot8(j) 

9 : Ci = Ci — t (Fij / Fipi vots (jj)Cpi vo t s (jj 

10: else 

11: Set pivot s(j) <— i 

12: Add to Z 

13: Break 

14: end if 

15: if Fi = then 

16: Add C 2 to Z 

17: end if 

18: end while 

19: end for 



Algorithm 3 Smith Normal Form Algorithm 

Input F :B Z 
pivots = 
for i = 1, . . . , n do 
while ^ do 

Set j lowest non-zero entry (smallest power of t) 
if pivots(j) ^ then 
A = diff(pivots(j),j) 

F t {Fj : i I ' Fjpi VO i s (j^Fpi VO t s ^j} 

else 

Set pivots(j) i 
for fc = j — 1, . . . , 1 do 
Y = diff(j,k) 
F* = F k -t Y {F k , i /F k , j )Fi 
end for 
end if 
end while 
end for 

All the pivots represent bars beginning at the grading of the corresponding generator and lasting the 
grading of the entry 

All zero columns represent infinite bars beginning at the grading of the corresponding generator 



there may be many mappings to get back to chain space. However, due to homogeneity, it is sufficient if 
we can compare two general elements of the row space and column space. Therefore in the above example, 
as we compute a basis, we would build an associated structure which given two boundaries or two cycles, 
returns their relative grading. 

Appendix B. Relative Persistent Homology Example 

Consider the 2-simplex, with each face introduced at a separate time. The extended homology of this is 
persistent homology of the persistent chains of this complex, with each simplex added to the torsion relations 
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in inverse order. Faces are introduced at times through 6, and added to the torsion relations at times 7 
through 13. 

The chain complex has the Unite presentation 

00 1 0"1) C2 1 (Toil cr 02, 012 > CTQ12 



c 



tcr i2 , £ 3 0i2 , t 5 0O2 , t 7 a 01 ,t 9 a 2 ,t 11 ai,t 13 ao 



The boundary operation defined on C is the normal boundary operation on simplicial complexes, defined 
by: 



<9(t = 
dai = 
dn 2 = 



d&oi = t 3 o"o — i 2 (Ti 

<?0O2 = ^ 4 0O — i 2 02 

9(T 12 = t 4 0i — t 3 a 2 

9CT012 = i 3 0Ol — £ 2 0O2 + ^012 



To compute the kernel module, we first compute its generators by doing Gaussian elimination on the matrix 
to the left, yielding the result to the right: 
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00 
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00 
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t 3 (T - t 2 (Tl 


001 







t 3 a - i 2 0i 


001 





i 4 cr - * 2 02 


002 







t 3 ai - t 2 a 2 


—tod + 002 





t 4 (Tl - t 3 (T 2 


012 










t 2 (Toi — to 2 +012 





1 — t 2 fT 2 + ta\2 


0012 







t 3 (T l - t 2 <7 02 + tU\2 


0012 





tooi2 





Pa 




— t(T012 





Po 


t 3 a 12 





Pi 




-t z a X 2 





Pi 


t 5 cr a2 





P2 




-i 5 O2 





P2 


t 7 a i 





P3 







£ 4 0O12 


P3 - tp2 + t- 


t 9 <r 2 





Pi 




-+ 9 2 





Pi 


t n (Tl 





P5 







i 7 0i 2 


P5 - tp4 


* 13 0o 





PS J 




V o 


i 9 0O2 


P6 - t 2 p4 



\ 



J 



This computation provides us with a module of generators with 6 elements, that injects into the gener- 
ators (To, . . . , 0012 to hit the elements <to, 0i, 02, i 2 0oi — too 2 + 012, £ 4 0oi2i i 7 0i2i i 9 0o2- We shall call the 
generators of the kernel module kq, . . . , kq with the images in the order stated here. They come in degrees 
0,1,2,5,10,12,13. 

To compute the relations for the kernel module, we setup another matrix for a Gaussian elimination; matrix 
on the left, result of the elimination on the right: 
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( 00 
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Ki 
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— ta a2 + 012 
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t 2 a i - ta Q2 + 012 
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— tooi2 





Po 
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-t 3 012 
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-i 3 0i2 





Pi 


-i 5 0O2 
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-£ 5 O2 
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-i 7 0oi 
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t 5 K 3 


P3 - tp2 + 


< 4 0O12 


K4 
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t 3 Po 


-i 9 2 





Pi 







i 9 K2 


Pi 


-t n (Tl 





P5 
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i 7 0i2 


K 5 












-t 4 pi 


-t 13 0o 
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t 13 n 


P6 


i 9 0O2 


K 6 
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K 6 


t 4 P 2 



t 2 pl 
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Reading off the kernel from the bottom of this reduced matrix gives us the complete presentation of the 
kernel module as: 

«o, • ■ • , «6 



ker d - 



t 5 K3, K4, t S K2,t 11 Ki, K§,t l3 Ko, Kq 



We can observe that by sheer luck, this presentation is already on Smith normal form, so no further work 
is needed. Furthermore, three of the relations - K4, k 5 , kq occur in the reduced matrix without any factor 
t, indicating they belong to barcodes of length - in fact, these are relations inherent in the generators at 
hand, but that were already in the torsion part of the module when they first materialize. This computation 
gives us the kernel module of the boundary map - the persistent cycle module. To compute homology, we 
now need to compute the cokernel of this module under the boundary map. First, we need to express the 
homology boundaries in the cycle basis. This can be done with yet another Gaussian reduction; result on 
the right: 



V 



00 


K 





01 


Ki 







K 2 





£ 2 0O1 — t&02 +012 


K3 





t 4 0~O12 


K4 





t 7 a 12 


K5 





t 9 <702 


Kq 











00 








01 








02 


-(tVo-tVl) 





001 


-(t 4 a a -t 2 a 2 ) 





002 


-(tVi-t 3 ^) 





012 


-(i 3 er i -£ 2 0O2 +ta 12 ) 





0012 


the cokernel, these images - 
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t 2 K2,t 4 Ki — t 3 K2,tK3 are included among the 



t 2 Ki, t 4 KQ 

relations, and a Smith normal form computed. The matrix for this Smith normal form is (0 replaced with • 
for clarity); again the computed Smith normal form is on the right: 



/t 3 
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t 5 



,13 



V 
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t 

v • 



13 



V 



From this normal form we can read off the final barcode for the homology module: 



Dim. 


Interval 





(0,13) 





(1.3) 





(2,4) 


1 


(5,6) 



